Quantum Trajectòries for Realistic Detection 
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Quantum trajectòries describe the stochastic evolution of an open quantum system conditioned 
on continuous monitoring of its output, such as by an ideal photodetector. Here we derive (non- 
Markovian) quantum trajectòries for realistic photodetection, including the effects of emciency, dead 
time, bandwidth, electrònic noise, and dark counts. We apply our theory to a realistic cavity QED 
scenario and investigate the impact of such detector imperfections on the conditional evolution of 
the system state. A practical theory of quantum trajectòries with realistic detection will be essential 
for experimental and technological applications of quantum feedback in many areas. 
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I. INTRODUCTION 



The limited utility of quantum measurement theory 
as axiomatized by von Neumann Q for describing prac- 
tical laboratory mcasurcmcnts has necessitated the de- 
velopment of more general measurement theories ^,|). 
In the past decade the application of such theories has 
become widespread in quantum òptics, in particular for 
describing continuous monitoring of the photoemission 
from radiatively damped open systems. They describe 
the evolution of the conditioned system state in terms of 
quantum jumps for direct detection and quantum 

diffusion ||,|7| for dyne detection. The stochastic evo- 
lution equation, termed a quantum trajectory, has also 
been applied in mesoscopic electrònics [|3|. 

Thus far, the main practical utility of quantum tra- 
jectory theory has been in improving the computational 
emciency of simulations used to compare models with 
experimental data. But it is now gaining increasing im- 
portance as the quantum generalization of Kalman fil- 
tering, which provides essential signal-processing meth- 
ods in classical estimation, communication, and control 
engineering. Quantum trajectory theory should in prin- 
ciple play the same pivotal role for emerging quantum 
analogs of these technologies || . Before this can hap- 
pen it is essential that the theory be extended to account 
for the imperfections of realistic measurement devices, as 
non-ideal detector dynamics can dramatically affect the 
proper inference from measured signals to the conditional 
quantum state of an observed system. 

In this paper we present the theory of quantum tra- 
jectòries for realistic photodetection. We model both 
photon countcrs and photoreceivers (for homodyne de- 
tection) and include the effects of efnciency, dead time, 
bandwidth, electrònic noise, and dark counts. The 
proper treatment of bandwidth limitations and electrònic 
noise are of particular significance as these imperfections 
are inevitable and predominant concerns in any practical 
context. They are of central importance in the current 



generation of experiments on quantum-limited measure- 
ment in atòmic and condensed matter flÏ3| systems. 

Our theory works by embedding the system within a 
supersystem that obeys a Markovian equation. If the set 
of (classical) detector states is S, then the supersystem 
is described by the set {p s : s £ S}. Here Tr[/9 S ] is the 
probability that the apparatus is in state s, and p s /Tr[p s ] 
is the system state given this event. 



II. THE SYSTEM 

In this paper we take the monitored system to be a two- 
level atom (TLA) , classically driven at Rabi frequency Í2 
and radiatively damped at rate T. The TLA obeys the 
unconditional master equation (ME) 

p = £p= -i(n/2) [a x ,p] + r (apj - \{J<x, p}) , (1) 

where a is the atòmic lowering operator, and a x = er + crí . 
Time arguments are not included unless they are neces- 
sary for the reader's understanding. In reality, it is diffi- 
cult to detect a significant fraction of an atom's fluores- 
cence. However, the ME ([!]) also describes, in a suitable 
regime |TJj], the damping of an atom through a cavity 
mode. This produces an easily detectable output beam. 
In this scenario, the effective decay rate T may be much 
larger than that of a bare atom, and we have this in mind 
when choosing F — 300Ms _1 for our simulations. 



III. PHOTON COUNTER 

An avalanche photodiode (APD) operating in Geiger 
mode produces a macroscopic current pulse in responsc 
to an incident photon. It consists of a p-n junction op- 
erated under a reverse bias greater than the breakdown 
voltage p5| . Under these conditions we can describe the 
diode by just three classical states (see Fig. [!]). The first 
(0) is a stable low-current state in which there are no 
charge carriers in the depletion region of the junction. 
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The transition from to the second state (1) takes place 
when an electron-hole pair is created in the depletion re- 
gion by an incident photon (with quantum efficiency rf) or 
by thermally initiated 'dark counts' occurring at a rate 
7dk- Further impact ionization leads to an avalanchc, 
until the current reaches some threshold value and a de- 
tection is registered, thus changing the state of the APD 
to 2. The transition from 1 to 2 has a random duration 
(we assume Poissonian) with mean 7 r _1 (the 'response 
time'). The avalanche is then arrested by the application 
of a negative-going voltage pulse that temporarily brings 
the bias voltage below the breakdown value fL7| |. This 
results in a fixed 'dead time', Tdd, during which the APD 
cannot detect photons, after which it is restored to state 
0. 

REALISTIC PHOTON COUNTER 
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FIG. 1. Realistic photon counting by an avalanche photo- 
diode. The quantum efficiency n is represented by the beam 
splitter (BS) . Single arrowheads within the realistic photode- 
tector indicate Poisson processes. For details, see text. 

Our aim is to derive the quantum trajectòries for the 
quantum system (the source of the light entering the 
APD) conditioned on the observation of an avalanche. 
For the TLA, we consider direct detection of fluores- 
cence. In this case the supersystem is described by the 
set (po,px,p2), where the tilde indicates that (for simplic- 
ity) we are using unnormalized system states, and the 
subscript indicates the associated detector states. The 
normalized conditioned TLA state is 



Pc = Pc/Tr^c 



Po + Pi + P2- 



(2) 



Our description of internal dynamics of an APD can be 
simply translated into rate equations for the discrete de- 
tector state s£ {0, 1, 2}, which in turn imply the follow- 
ing stochastic generalization of the ME ([]]): 

dp Q = dt { \C - 7 dk - rpJ - AT] po + N{t*)~P2} , (3) 
Pi + V^Jpo + 7dk/5o| , (4) 
P2+Npi}- (5) 



dpi = dt | C — 7 r — M 

Hc-Àí(t*) 



dp2 = dt ■ 



Here Jp$ denotes cpocrL We use Aí for the number 
of detections counted, so that dM (t) = Mdt is a point 
process equal to 1 in the infinitesimal interval when an 
avalanche is first observed and otherwise. The delayed 
process, M(t*) = M(t — Tdd), is used to return the de- 



tector to state 0. The statistics of dM are defined by its 
expectation value E[dA/] = 7 r díTr[ ( 5i]/Tr[ j 5 c ]. 

The detector imperfections lead to substantial changes 
in the conditional dynamics of the TLA, as compared to 
ideal quantum trajectòries. Representative features can 
be seen in Fig. || (A) and (B). Plot (A) shows a typi- 
cal portion of a trajectory for z c , while plot (B) shows 
the same, and y c , over a shorter time around t sa 4.9, 
when an avalanche is registered. Unlike the case of ideal 
detection, the corresponding "quantum jump" does not 
take z c —>■ — l,y c ~ * 0, and the amplitude of subsequent 
oscillations in z c ,y c are less than 1. The jumps in the 
conditioned quantum state caused by the detection of 
avalanches are attenuated because of the finite detector 
response time in combination with the continuous Rabi 
oscillation, which evolves the TLA away from the ground 
state for a random and unknown time (with mean 7, ) 
between the "actual" spontaneous emission event and the 
registration of the photocurrent avalanche. During the 
APD dead time the effective efficiency is zero, and as a 
result the TLA's conditional state regresses towards the 
steady state of the unconditional ME (1). Even after the 
detector becomes ready again (by resetting to state 0) 
the Rabi oscillations in z c , y c decay because of the APD's 
non-unit efficiency and finite bandwidth. 

These imperfections cause the stationary ensemble- 
averaged conditional purity p = limt^oo E {Tr[p^(í)]} to 
be substantially less than one for large íl. For small 
Í2 however, even the unconditional (without measure- 
ment) stationary purity p u of the TLA approaches unity. 
It is thus useful to define a sealed purity e [0, 1] that 
measures how much improvement measurement gives: 
Sealed p — (p — p u )/(l — p u ). For the typical parameter 
vàlues used in Fig. 0, the Sealed p ~ 0.052. 




Time (Units of r _1 ) 

FIG. 2. In plot (A), z c = Tr[(2CT f CT - l)p c ] is shown for 
a typical trajectory of duration 20r~ 1 . In plot (B), z c and 
y c — — iTr[((j — a^)p c ] (dash-dot) are shown for the same tra- 
jectory near the time of the second APD avalanche. The real- 
istic parameters used for this photodetection trajectory were 
n = 80%, and (in units of the TLA decay rate F = 300Ms _1 ) 
7r = 7, r dd = 2, 7<jk = 5 x 10~ 6 and fi = 10. 
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IV. PHOTORECEIVER 

When the incident photon flux is high, as in homo- 
dyne detection, a p-i-n photodiode connected to a tran- 
simpedance amplifier (see Fig. ||) is an appropriate pho- 
toreceiver 17 . When a photon strikes the depletion rc- 



gion of the p-i-n junction, an electron-hole pair is pro- 
duced, with probability equal to the quantum efficiency 
77. The charge carriers drift under the influence of the 
below-breakdown reverse bias, and the resultant current 
I is fed into an operational amplifier (op-amp) set up as a 
transimpedance amplifier. This has a low effective input 
impedance, so that the diode acts as a current source, 
and I is converted into a voltage drop V across the feed- 
back resistor, R. The capacitor C, in paral·lel with R, 
represents the total capacitance from the output of the 
op-amp back to its input, including capacitance added 
deliberately for the smoothing of noise and oscillations. 
If no electrònic noise were present, the output voltage of 
the photoreceiver would be a filtered version of the input 
signal given, in the frequency domain, by 



V(uj) = -IR/(l + iuRC). 



(0) 




REALISTIC PHOTORECEIVER 
FIG. 3. Homodyne detection by a photoreceiver. The out- 
put field of the TLA is combined with a LO before being 
detected by a realistic photoreceiver consisting of a p-i-n pho- 
todiode (of quantum efficiency rf) that produces the photocur- 
rent 7, and an ideal op-amp with feedback resistor R, and 
capacitance, C (see text). The output voltage V is what is 
measured in the laboratory. 

It should be noted that if this were the case (that is, if 
there were no noise) then the input / could be perfectly 
reconstructed from the filtered signal V. Thus the re- 
sultant quantum trajectòries would be no different from 
those of a photoreceiver with infinite bandwidth. Ev- 
crything of interest results therefore from the presence 
of excess noise. We include only the Johnson noise Vj 
from the feedback resistor, which has a flat spectrum 
Sj = Ak^TR. This simplification (ncglecting contribu- 
tions from voltage noise of the operational amplifier) can 
be justified for practical receivers with R ~ lOkíl. 

The output voltage V from the photoreceiver is given 
by the sum of the filtered signal and the Johnson noise 



V = V + Vj. 



(7) 



Our aim is to find the quantum trajectory for the sys- 
tem, conditioncd on continuously monitoring V. Sincc 



the voltage V, which describes the detector state, is a 
continuous variable, in this case S = 1R, the real line, and 
the supersystem can be described by an operator func- 
tion p(V). Finding the stochastic equation of motion for 
p(V) is quite involved. 

We begin by taking the output current I of the photo- 
diode to be that from a perfect (apart from its efficiency 
77) unbalanced homodyne detection of the fluorescence of 
the TLA. For a LO tuned to the atòmic transition fre- 
quency u>o, of power P, and phase (f>, the current is [|6jjl9|] 



I = e^fpjh7 Q WiVfie-^a + e'V) + y/rjÇ(t) , (8) 

where we have ignored the D.C. component due to the 
LO power. Here £(í) is the Gaussian white noise jÏ8j aris- 
ing from the Poissonian statistics of the LO and e is the 
electron charge. The evolution of the TLA conditioned 
on I is given, in terms of the noise by the following 
stochastic master equation |19|| 



(9) 



dpi = dt {£ + y fijrÇ(t)H[e- i +a] } p, 



where H[A]p = Ap + pA^ - Tr[Ap + pA^]p. 

Now Eq. (^|) is equivalent to the stochastic equation 



/ + V/R + C(dV/dt) = 0. 



(10) 



Since the voltage V is not directly measured, we must 
consider a distribution P(V) for it. Assuming that 
C > 0, and, for the moment, that I is known, Eq. ( p~0| ) 
can be converted to an Itò [Q stochastic Fokker-Planck 
equation for P(V) conditioned on the photocurrent: 



dPi(V) = 



í d V + IR ^ Prje 2 
\dV RC 



d 2 



2Hlü C 2 dV 2 



P(V)dt. (11) 



Here we are using the convention that subscripts indicate 
that the increment is conditioned on that result. That 
is, for example, P^V) = P(V\I). 

Next we need to determine the effect of the measure- 
ment of V on P(V). This can be calculated by using 
Bayes' conditional probability theorem 



P V (V)=P V (V)P(V)/P(V). 



(12) 



Remembering that the Johnson noise is white, it follows 
from Eq. ([?]) that Pv(V) is a Gaussian with mean V and 
variance Ak-oTR/ dt. From this we find that 



P(V) = J dVP v (V)P(V). 



(13) 



is a Gaussian of mean (V) and variance Ak-gTR/ dt. It 
follows that we can writc 



V = (V) + ^4,k B TRdWj(t)/dt, 



(14) 
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where dWj(t)/dt is another Gaussian white noise source, 
independent of £(t). Substitution of -Py(V) and P(V) 
into Eq. ( |ï^ ) yields the effect of the V-measurement: 

dP v (V) = {V- (V)) P{V)dWj{t)/^jAk B TR. (15) 

Now, to see how V conditions the TLA, we form the 
quantity p(V) = pP(V), where p is here independent 
of P(V) because we are imagining I to be known at all 
times. The time evolution of p(V), given that V and I 
are known, is found from 

p(V) + d Piy (V) = (p + d Pl ) 

x [P(V)+dP!(V)+dP v (V)], (16) 

with the use of Eq. (jï|), Eq. |ï|) and Eq. (§). 

Finally, in reality, V is known but / is not. There- 
fore we should average over the vacuum noise £(í), but 
keep the Johnson noise dWj jdt. We define a dimension- 
less voltage v = V y/C/Ak^T , a rate 7 = 1/RC and a 
dimensionless noise power N — Ak^TfiujQ/rjRPe 2 . This 
last expression is the ratio of the low-frequency power 
in V from the Johnson noise to that from the vacuum 
noise. We then obtain the following stochastic nonlincar 
superoperator Fokker-Planck equation for p(v): 

( 7 <9 2 d \ 
dp v (v) = dt{£+—-s+ 7 -v) p(v) 

+ dt^ [e-**vp{v) + e^p(v)aï] 

+ ^dWj(t)(v~(v))p(v). (17) 

The dependence on V may be explicated by substituting 
dtj (y/C/4k B TV - (v)) for ^fdWj(t) [see Eq. 0)]. In 

the above, (v) = J dvTr[p(v)]v. The normalized condi- 
tioned TLA state is p c = J p(v)dv. 

A Typical trajectory for realistic homodyne x ((f> = 0) 
detection is shown in Fig. |^ (A) . The main difference from 
the case of perfect detection is the reduced amplitude of 
variation in x c (dotted) and z c (sòlid) . This is due to the 
effective bandwidth of the photoreceiver, which affects z c 
more because of its faster dynamics. Plot (B) shows the 
photoreceiver output voltage V that is used in Eq. ( |Ï7| ) to 
condition the TLA state. It is seen that V is correlated 
with expected. 

See attached flle TrajHom.jpg 

FIG. 4. Plot (A) shows x c = Tr[(cr + <T t )p ] (dotted) and 
z c (sòlid) for a realistic homodyne x trajectory. Plot (B) is 
the dimensionless output voltage from the photoreceiver. The 
photoreceiver parameters were N = 0.1, 77 = 98%, 7 = 1.5F 
System parameters were as for Fig. ^| Plot C gives the Scaled 
p as a function of the driving, Q, for homodyne x (dotted) and 
y (dash-dot) detection. 



Plotted in (C) is the scaled purity as a function of the 
driving strength for both homodyne x and y (<f> — 7r/2) 
detection. As íl increases, homodyne y detection be- 
comes increasingly worse than x detection at follow- 
ing the evolution of the TLA. This is due to the finite 
bandwidth of the photoreceiver in combination with the 
conditional homodyne dynamics in the Í2 3> T limit 
0. For homodyne x measurement the x quadrature, 
which changes sign fairly infrcqucntly, dominates the 
TLA state. The slow (r) dynamics allow the detector 
to track of the state reasonably well. In contrast, homo- 
dyne y detection produces a conditional state dominated 
by fast (Í2) Rabi cycling, which is poorly followed. 

V. CONCLUSIONS 

In conclusion, we have presented a theory of quantum 
trajectòries for systems conditioned on realistic photode- 
tection. The equations are tractable, as we have demon- 
strated by numerical simulations, and allow us to quan- 
tify the degree and manner by which imperfections such 
as a finite bandwidth modify the conditioning of quan- 
tum states by measurement in concrete experimental sce- 
narios. Realistic quantum trajectory models will be of 
paramount importance in the field of real-time quantum 
feedback control pc| ]. The techniques we introduce here 
may also prové essential in describing other realistic mea- 
surements, such as in condensed matter systems |||. 
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